EJPE 


European J of Physics Education Vol. 3 Issue 3 2012 


Hu 


Deriving the Work Done by an Inverse Square Force in Non-Calculus-Based Introductory 

Physics Courses 


Ben Yu-Kuang Hu 

Department of Physics, University of Akron 
Akron, OH 44325-4001, USA 
byhu@uakron.edu 


Abstract 


I describe a method of evaluating the integral of V 2 with respect to r that uses only algebra and the concept of area 

underneath a curve, and which does not formally employ any calculus. This is useful for algebra-based introductory 
physics classes (where the use of calculus is forbidden) to derive the work done by the force of one point charge or 

mass on another. I extend the method to evaluate any integral of the form J'r" dr for any integer n ^ — 1. 
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Introduction 

In the United States, a large proportion [approximately 40% in 2008 according to Mulvey and 
Nicholson (2011)] of university students who enroll in introductory physics courses take the 
algebra-based version, in which formal calculus is not employed. The popularity of the algebra- 
based courses is in part due to typical medical school requirements and the Medical College 
Admissions Tests syllabus (American Association of Medical Colleges (2008)), for which 
algebra-based version of the introductory courses satisfies the physics requirement. 

In algebra-based introductory physics courses, the result for the work done by an inverse 
square force, such as the Coulomb or gravitational force between point particles, is usually given 
without proof, since the rigorous derivation of the result requires knowledge of integral calculus. 
Here, I describe a way that the work can obtained without any formal calculus, based on work 
being the area under the force vs. displacement curve. 

The derivation for area under a l/ r 2 curve 

The concept of work in algebra-based introductory physics courses is typically introduced for the 
special case of a constant force, in which case the work done is W = Fd, where F is the 
component of the constant force in the direction of the displacement and d is the distance 
travelled. This corresponds to the area enclosed in the rectangle created by the constant force vs. 
displacement curve. The work done by a variable force is later explained as the area underneath 
the force vs. displacement curve, typically with the aid of a figure such as shown in Fig. 1, where 
the area underneath the curve has been divided into many small rectangular blocks. In the case of 
a Hooke’s law force which is proportional to the displacement, F = kx, the area underneath the 
curve between x, and xp is IF = fk(x 2 - xf ). This can be obtained without calculus because it is 
the area of a right-angled trapezoid, which is the sum of the area of a rectangle and a right-angled 
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triangle. However, in the case work done by forces that vary as the inverse-square of the 
distance, F = k/r 2 , the work done in moving a particle from r, to rf, 


W = k£dr\=k 


\ 


'f) 


( 1 ) 


is usually simply stated as a result whose derivation requires knowledge of mathematics that is 
beyond the scope of a non-calculus based course. Consequently, students often do not remember 
this result. In my experience, when students are asked in tests for the work done by a force 
F = k/r 2 in moving an object from a distance R to R+d, instead of the correct answer 


W = k\-- 


1 


R R + d 


, a significant fraction of the students recall the constant force result and 


answer W = Fd, where they typically use the initial force k/R 2 for F. This error is less common 


when a similar question is asked for a Hooke's law force, in my opinion, because students are 
explicitly shown that the work done for a spring is that of the area underneath the curve, whereas 
for the inverse-square forces, they are simply told the result as a fact. If this is the case, the 
solution is to derive the result for them in the inverse-square force case. But how do we do this 
without the use of integral calculus, and in a simple and straightforward way? 



Figure 1: Area underneath a curve F from r, to rf, divided into N subintervals, 
where the limits of sub-intervals are r t = r 0 , r lt r 2 , ... , r N _ It r N = ry. The area Aj 
under the curve in the sub-interval from rjand rj_i is equal to (rj— rj_i) <F>, where 
<F> is the average of the function F over that interval. 
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The integral ( dr — is the area underneath the — curve from r, to rf. Divide the interval 

Jl; y- y 1 

q to rf into N equally spaced sub-intervals, where N is a large number and the j-th sub-interval is 
from rj-i to rj (with ro = r, and r\ = rf), as shown in Fig. 1. The area under the curve in the j-th 
sub-interval from rj_i to rj is 


A j = (r - r j-i)(jjj ( 2 ) 

where / —-\ is the average of —- over the sub-interval. Since the number of intervals N is large, 

V / r 

the difference between rj_i and rj is very small and we can substitute any r from the interval ryi 
and r, into the r in /-^\and obtain an accurate result. Let us choose to substitute (with the 


benefit of hindsight!) the geometric mean of the lower and upper limits of the sub-interval, 

‘ 1 \ 1 


■( oo -.) I/2 


. Thus, ( — 


00-i 


which gives 


r ■ 

A i = ~ 


- r 


j -1 


0 0-1 


1 1 


'j -1 


0 


(3) 


The area underneath the — curve from q to rf is the sum of the areas underneath the 

r 

curves in the sub-intervals r 0 (= q) to q, ri to r 2 , r 2 to q, and so on until r N -i to r N (= q-), which 
gives 


.. 1 -V -V 1 1 

1 ./ = 1 7 = 1 0-1 0 


i a n 


j \ 

i 


0 r \ 

l 1 


(\ l\ 


r x 0 / 


+ ■■• + 


/ 1 1 \ 


r N-2 V N -1 / 


(l 1 \ 


\ r N -1 r N ) 


0 r N 


r. r, 


(4) 


where the second-to-last equality results from second term in the first parenthesis cancelling the 
first term in the second parenthesis, the second term in the second parenthesis cancelling the first 
tenn in the third parenthesis, and so on. 


Generalization to Jr" dr for integers n * 1 

The above result can be generalized to integrate r 11 for any integer n, except n = -1. As in the 

case of integrating \ , one divides the integration interval from q to rf into sub-intervals given 

y~ 
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by r 0 = r,; n; r 2 ; • • ■ ; r N -i; r N = rf, where N is a large number, so that the sub-intervals are very 
small. I deal with the cases —and r 11 for n > 1 separately. (The r° case is trivial because r° = 1.) 


Evaluating Jdr for integers n > 1 


The area under the curve —- in the sub-interval from r,_i to r, is 

r „+i j j 

a j = (0 - r j- ( 5 ) 

where (\/r n+l ^j is the average over the interval. Because the sub-interval is so small, we can 

replace l/r 11+1 with any function of q.i and r J? which gives l/r 11+1 when both r,.| and q are replaced 
by r. Again, with the benefit of hindsight, we choose 



r /~‘ + 0”~ 2 Q-i + r j~ 3r ji + • • • + r j r 'j-\ + r 'j-\ 

r]rU 


( 6 ) 


(Note that there are n terms in the numerator on the right hand side of Eq. (6), so when q and q_i 
are replaced by r, the numerator becomes nr 11 ' 1 .) With this choice, 

' i \ (r -'}-■) ( r r' +r r 2r j-i +r r' r l 


4 -(0 -rj_M 




;) 


n+\ 


rJr'U 


The numerator in Eq. (7) “telescopes”; i.e., 


(V) 


/ \/ n -1 , n-2 , n -3 

(Tj-rj-iW +r J r i-i + r j r j 


M + ... + r jr ;:{+ rj _ 


n -1 | n 

M ) =r j 


, n -1 , n-2 2 , , n -1 

+ V J 0-1 + 0 0-1 +••• + 00-1 


"(o 


n -1 


n-2 2 


J 0-1 + o- ~0-i + • • • + 00-1 + 0-1) = 0" - 0-1 


( 8 ) 


and Eq. (7) becomes 
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4 =- 

n r 


( n n \ 

0 “O-i/ 


j 0-1 


1 1 
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1 / 


( 9 ) 


Therefore, 
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Evaluating Jr' 1 dr for integers n > 1 

The area under the curve r " in the sub-interval from rj_i to rj is A. = (ir - r._, )^r"^. With the 
benefit of hindsight, we replace 

f")-Wf; + ^_, + ^, + ... + ^' + ^) (12) 

(Note that there are n+1 terms in the parentheses of the above equation, which gives (n+l)r n 
when rj and r,.| are replaced by r.) Using the “telescoping” result, Eq. (8), gives 


4 - (r - oW ') - f+'TV. + 

' ' 77 + 1 


'fN, + ... + r . 


rd+rj) -(13) 


77 + 1 


and therefore 

r f N 


C i n V'' a 1 V'' n + 1 n+1 1 \( ii + l n+1 1 . ( n + 1 n + 1 1 . . ( n + 1 n + 1 \ . ( n+1 ,i+l M 

fdrr = \ A = -- \ r - r = - -[{r { -r 0 + [r 2 -r x + ■ ■ ■ + (r v _, - r N _ 2 + [r N - r N _ x J 

J r . fa n + lfa n +1 


--i-k'-cl-Wk 1 -<"] ■ 

77 + 1 77 + 1 


(14) 


Conclusion 

In conclusion, for the special but important case of the work of an inverse-square force F = k/r 2 , 
one can demonstrate that the work done by the force in moving a particle from r, to rf is 
k[y - y- jwithout the formal use of integral calculus. Many students in non-calculus-based 

introductory physics classes have a fixation with the constant force equation for work, W = Fd, 
and will use it even when the force is not constant. Therefore, an explicit demonstration of the 
result for the inverse-square force case should be helpful in preventing these students from 
making that common mistake. The method described is also extended to evaluate the integral of 
r n with respect to r for any integer 77 ^ 1 . 
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